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We consider the quantum and classical Liouville dynamics of a non-integrable model of two 
coupled spins. Initially localised quantum states spread exponentially to the system dimension 
when the classical dynamics are chaotic. The long-time behaviour of the quantum probability dis- 
tributions and, in particular, the parameter-dependent rates of relaxation to the equilibrium state 
are surprisingly well approximated by the classical Liouville mechanics even for small quantum 
numbers. As the accessible classical phase space becomes predominantly chaotic, the classical and 
quantum probability equilibrium configurations approach the microcanonical distribution, although 
the quantum equilibrium distributions exhibit characteristic 'minimum' fluctuations away from the 
microcanonical state. The magnitudes of the quantum-classical differences arising from the equilib- 
rium quantum fluctuations are studied for both pure and mixed (dynamically entangled) quantum 
■ states. In both cases the standard deviation of these fluctuations decreases as {Ti/ JY^^ , where J 

is a measure of the system size. In conclusion, under a variety of conditions the differences between 
quantum and classical Liouville mechanics are shown to become vanishingly small in the classical 
limit {J /Ti — > oo) of a non-dissipative model endowed with only a few degrees of freedom. 
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lO ■ I- INTRODUCTION 

o 

' The study of chaos in quantum dynamics has led to differing views on the conditions required for demonstrating 
quantum-classical correspondence Moreover, the criteria by which this correspondence should be measured 

have also been a subject of some controversy While much of the earlier work on this topic is concerned with 

characterizing the degree of correspondence between quantum expectation values and classical dynamical variables 
[^-H, the more recent approach is to focus on differences between the properties of quantum states and associated 
classical phase space densities evolved according to Liouville's equation [p]j9|-p^]. 
^ , Several authors have examined quantum-classical correspondence by considering the effects of interactions with a 
' stochastic environment [^-16|, a process sometimes called decoherence. While this process may improve the degree 
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^ : of quantum-classical correspondence for fixed quantum numbers, it has been further suggested that the limit of large 
O^' quantum numbers is inadequate for correspondence, and that decoherence must be taken into account to generate 
classical appearances from quantum theory; this view has been argued to apply even in the case of macroscopic bodies 
, , that are described initially by well-localised states, provided their classical motion is chaotic ||,|l^,|I^ . In this paper 
r^J ' we examine how the degree of correspondence with Liouville dynamics scales specifically in the limit of large quantum 
^ . numbers. This "classical limit" is distinct from a "thermodynamic limit" , that is, a limit involving many quantum 
numbers. 

The degree of Liouville correspondence has been characterized previously by studying the differences between the 
means and variances of the dynamical variables ||l|,p|,pd]- |r3| , p^ . This involves a comparison of quantum expectation 
values and classical ensemble averages. However, these low-order moments give only crude information about the 
differences between the quantum and classical states. Specifically, the quantum state may exhibit coarse structure 
which differs significantly from the classical state although the means and variances (for some simple observabes) 
are nearly the same for the quantum and classical states. Moreover, much of the previous work was concerned with 
correspondence at early times, or more precisely, in the Ehrenfest regime when the states are narrow compared to 
system dimensions p^ , p^ . 

Another approach is to identify quantum-classical differences with differences between the Wigner quasi-distribution 
and the classical phase space density This approach is objectionable because the Wigner quasi-distribution may 
take on negative values and therefore may not be interpreted as a "classically observable" phase space distribution. 
It is possible to consider instead smoothed quantum phase space distributions, but in this case the residual quantum- 
classical differences still do not have clear experimental significance. 

In this paper, we characterize the degree of quantum-classical correspondence by comparing quantum probability 
distributions for dynamical variables with the corresponding classical marginal distributions for these dynamical 
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variables. These are well-defined classical observables that describe the distribution of outcomes upon measurement 
of the given dynamical variable. We are interested in the differences that arise on a fine scale and therefore characterize 
the typical quantum-classical deviations that arise in bins of width h. 

The dynamics are generated by a model of interacting spins that we have studied previously in Ref. ||l3|] . The 
Hilbcrt space is finite dimensional so no artificial truncation of the state is required. The quantum time-evolution 
is unitary and the classical motion is volume-preserving (symplectic). In the case of classically chaotic motion, we 
follow initially localised states until they have evolved well beyond the relaxation time-scale of the classical density. 
Throughout the paper we emphasize that the quantum signatures of chaos that appear in the quantum distributions 
are the same as those that appear in the marginal classical distributions. In particular, the quantum relaxation rates 
can be accurately estimated from the Liouville dynamics of an approximately matching initial phase space density. 
This purely classical approximation is surprisingly accurate even for small quantum numbers, but may be most useful 
for the theoretical description of mesoscopic systems since the purely classical calculations do not scale with the 
quantum numbers. 

The quantum and classical probability distributions remain close even after the states have spread to the system 
dimension. Specifically, in mixed regimes, the quantum distributions exhibit an equilibrium shape that reflects the 
details of the classical KAM surfaces. When the classical manifold is predominantly chaotic, the quantum and classical 
states relax close to the microcanonical state. However, in both of these chaotic regimes the equilibrium quantum 
distributions exhibit characteristic fluctuations away from the classical ones. We demonstrate that the standard 
deviation of these quantum-classical differences becomes vanishingly small in the classical limit, /h — *■ oo, where J7 
is a characteristic system dimension. 

This paper is organised as follows. In section II we describe the quantum and classical models for our system. 
In section III we describe the initial quantum states and corresponding classical densities. In section IV we give an 
overview of the dynamics of the probability distributions in the different classical regimes. In section V we examine 
the quantum and classical relaxation rates using the Shannon entropy. In section VI we give an overview of the 
time-development of the small differences between the quantum and classical probability distributions. In section VII 
we show that the relative quantum-classical differences decrease as an inverse power of the Hilbert space dimension. 
In section VIII we provide a brief summary and explain how our results inform current discussions of the necessary 
and sufficient conditions for the emergence of classical properties from quantum mechanics. 

II. THE MODEL 

We consider the quantum and classical dynamics generated by a non-integrable model of two interacting spins, 

oo 

H = a{S, + L,) + cS^L^ ^ 5{t - n) (1) 

n— — oc 

where S = {Sx, Sy, Sz) and L = (ij,, Ly, Lz). The first two terms correspond to simple rotation of both spins about 
the z-axis through an angle a with range 27r radians. The sum over coupling terms describes an infinite sequence of 
(5-function interactions at times t = n for integer n. Each interaction term produces an impulsive rotation of each 
spin about the x-axis by an angle proportional to the a;-component of the other spin. 

A. The Quantum Dynamics 

To obtain the quantum dynamics we interpret the Cartesian components of the spins as operators satisfying the 
usual angular momentum commutation relations, 

[Si, Sj] — ieijkSk 
[Li, Lj] = ieijkLj; 
[Ji, Jj] — ieijkJk- 

In the above we have set Ti ^ 1 and introduced the total angular momentum vector J = S + L. 

The Hamiltonian (0) possesses kinematic constants of the motion, [S^, H] = and [L^, H] = 0. Thus the total state 
vector [ip) can be represented in an invariant Hilbert space H = Hs ^ Ti-i, with dimension N = (2s -I- 1) x {21 + 1), 
that is spanned by the orthonormal vectors 

|s, I, rus, mi) = |s, ms) (g) 1^, mi) (2) 
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with nis G {s, s — 1, . . . , —s} and m; G {/, / — 1, . . . , —I}. 

It should be noted that the components of the total angular momentum are not conserved [Ji^H] ^ 0. The z- 
component is subject to the selection rule, A = {±2, 0}, and consequently the full Hilbert space can be decomposed 
into two invariant subspaces. 

The periodic sequence of interactions introduced by the (5-function produces a quantum mapping. The time- 
evolution for a single iteration, from just before a kick to just before the next, is produced by the unitary transfor- 
mation, 

\i>{n + l))=F\i,{n)}, (3) 

where F is the single-step Floquet operator, 

F = exp [-ia{S^ + L^)] exp [-icS^ (g) L^] . (4) 
The quantum dynamics are thus specified by two parameters, a and c, and two quantum numbers, s and I. 



B. Classical Map 



For the Hamiltonian ([l|) the corresponding classical equations of motion are obtained by interpreting the angular 
momentum components as dynamical variables satisfying. 



{Si, Sj} - 
{Li, Lj} 



^ijkLk 
^ijkJk, 



with {•,•} denoting the Poisson bracket. The periodic J-function in the coupling term can be used to reduce the 
time-evolution to a stroboscopic mapping at times t — n, for integer n. 



cn+l 
cn+l 



f n+1 



'n+1 _ 



cos(a) - 5^ cos(7rL;') - SI' sinijrV^) sin(a), 

S^ cos{jrL2) - S^ sin(7rL^)j cos(a) -I- S!^ sin(a), 
^:cos(7rL::) + 5,"sin(7rZ::), 

cos(a) - [l; cosi-fS^) - sin(75:)] sin(a), 
'l; cos(74") - ^" sin(74")] cos(a) + sin(a), 
L>os(7^:) + Z,^sin(7^:). 



(5) 



Here L = L/|L| , S = S/|S|, and we have introduced the parameters 7 — c|S| and r = |L|/|S|. The mapping equations 
describe the time-evolution of from just before one kick to just before the next. 

Since the magnitudes of both spins are conserved, {S^,H} = {L^,H} = 0, the stroboscopic motion is actually 
confined to the four-dimensional manifold V = x S"^, which corresponds to the surfaces of two spheres. This is 
manifest when the mapping (j^) is expressed in terms of the four canonical coordinates x = (Sz, Lz, (f>i), where 
(ps = ta.n.{Sy/ Sx) and (pi = taii{Ly/Lx). We will refer to the mapping ^ in canonical form using the shorthand 
notation x"+^ F(x"). It is also useful to introduce a complete set of spherical coordinates 9 = {9s,4's,Si,<f>i) where 
Os = cos-i(S',/|S|) and Oi = cos~i(i./|L|). 

As in the quantum model, the components of total angular momentum are not constants of the motion { J^, H} ^ 0. 
On the other hand, the quantum selection rule AJ^ — {±2,0} has no classical analogue. 

The mapping on the reduced surface V enjoys a rather large parameter space. The dynamics are determined 
from three independent dimensionless parameters (a, 7, and r > 1), where 7 = c|S| is a dimensionless coupling 
strength and r = |L|/|S| corresponds to the ratio of the magnitudes of the two spins. The dependence of the classical 
behaviour on these parameters is described in Ref. [[l3[. 
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C. The Liouville Dynamics 



We are interested in comparing the quantum dynamics generated by (|^) with the corresponding LiouviUe dynamics 
of a classical distribution. The time-evolution of a Liouville density is generated by the partial differential equation, 

g-^ = -{Pc,H}, (6) 

where H stands for the Hamiltonian (Q) and x = {Sz, (fis, Lz,4>i)- 
The solution to (^) can be expressed in the compact form, 

Pe(x,t)= / dA*(y)<5(x-x(i,y))p,(y,0), (7) 
J-P 

with measure, 

d^y) = dSzd(j)sdL,d(j)u (8) 

and where the time-dependent functions x(<, y) € 7^ are solutions to the equations of motion (|^) with initial conditions 
y G V. This solution expresses the fact that Liouville's equation (||) describes the dynamics of a classical density of 
points evolving in phase space under the Hamiltonian flow. We exploit this fact to numerically solve (|^) by randomly 
generating initial conditions consistent with an initial phase space distribution /3c(x, 0) and then time-evolving each 
of these initial conditions using the equations of motion (W) . 



D. Correspondence Between Quantum and Classical Models 

For a quantum system specified by the four numbers {a, c, s, 1} we determine the corresponding classical parameters 
{a, 7, r} by first defining the classical magnitudes in terms of the quantum magnitudes, 

|S| = v^^Ti) 

|L| = ViifTT), (9) 

where the quantities on the left hand side are the lengths of the classical spins and those on the right are the quantum 
numbers. If we set the Hamiltonian coefficients a and c numerically equal for both models, then the remaining two 
dimensionless classical parameters are determined, 

V + 

7 = c^/s{s + l). (10) 

We are interested in extrapolating the behaviour of the quantum dynamics in the limit s — > oo and I oo. This is 
accomplished by studying sequences of quantum models with increasing s and I chosen so that r and 7 are held fixed. 
Since s and I are restricted to integer (or half-integer) values, the corresponding classical r will actually vary slightly 
for each member of this sequence (although 7 can be matched exactly by varying the quantum parameter c slightly). 
In the limit s — > 00 and / — > 00 this variation becomes increasingly small since r = + l)/s(s -I- 1) l/s. We 
have examined the effect of the small variations in the value of r on the classical behaviour and found this variation 
to have negligible effect. 

III. INITIAL STATES 

We consider initial quantum states which are pure and separable, 

\m) = \M0))^\M0))- (11) 

The initial state of each subsystem is a directed angular momentum state, 
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\e,<j,) = R(i\e,4>)\j,j), 



(12) 



where j in this section refers to either / or s. This is a localised state, i.e. one of maximum polarization in the direction 
(0, 0), with expectation values of the spin components confined to the surface of a 2-sphere, 



{e,cj)\j,\e,4,)^3cose 



\Jx i lJi 



= je^*"^ sin( 



(13) 



The states (g^ are the SU(2) coherent states, which, like their counterparts in the Euclidean phase space, are 
minimum uncertainty states pO[; the normalized variance of the quadratic operator. 



^|J2|0,, 



JJ|^,, 



is minimised for given j and vanishes in the limit j oo. The coherent states directed along the z-axis, 
|j, — j), saturate the inequality of the Heisenberg uncertainty relation. 



{Jl){jy)> 



{Jz 



(14) 



and 



(15) 



although this inequality is not saturated for coherent states polarized in other directions. 

We would like to construct a classical Liouville density on the 2-sphere with marginal distributions that match 
the quantum probability distributions. But we have shown previously that this is impossible for the SU(2) coherent 
states |l^. Thus from the outset it is clear that any choice of initial classical state will exhibit residual discrepancy 
in matching some of the initial quantum moments. 

We have examined the correspondence properties of several different classical distributions. These included the 
vector model distribution described in the Appendix of and the Gaussian distribution used by Fox and Elston in 
correspondence studies of the kicked top 110] . We selected the density. 



sin 9d9d(p ~ C cxp 



2sin2(i)) 



sin 9d9d(p 



(16) 



with C = [27rcr^ (l — exp(— 2cr~^))] , instead of those previously considered because it is periodic under 2n rotation. 
The classical density (^6|) has a maximum along the -|-2;-axis, corresponding to the coherent state |j, j). An initial 
state directed along {9o, (j^o) is produced by a rigid body rotation of (16) by an angle 6o about the y-axis followed by 
rotation through an angle 4>o about the z-axis. The variance is a free parameter of the distribution. Although 
cannot be chosen so that all low order moments are satisfied, the choice = 2|J|, where |Jp = j{j + 1), produces 
a reasonable compromise, as discussed in [O. 



IV. DYNAMICAL BEHAVIOUR OF PROBABILITY DISTRIBUTIONS 

In the case of a mixing classical system, initial densities with non-zero measure are expected to spread in an 
increasingly uniform manner throughout the accessible phase space. The term uniform is meant to apply specifcally 
in a coarse-grained sense. For some simple maps, such as the baker's map, it is possible to show that this rate of 
relaxation to the equilibrium configuration occurs exponentially with time [pi| . 

The spin map we consider (^ is not mixing on the accessible classical manifold V, but has mixed dynamics: 
depending on the system parameters, the surface V can generally be decomposed into regions of regular dynamics 
and a connected region of chaotic dynamics [^3|. In parameter regimes that are predominantly chaotic, we expect 
behaviour on V that approximates that of a mixing system. In particular, initially localised Liouville densities should 
relax close towards the microcanonical measure at an exponential rate, on average. In this section we demonstrate 
that these signatures of chaos are exhibited also by the quantum dynamics. Most striking is the degree of similarity 
between the quantum and classical behaviours even in regimes with classically mixed dynamics. 

We are interested in the behaviour of quantum probability distributions that are associated with measurements 
of classical dynamical variables. The quantum probability distribution associated with the classical observable is 
given by, 



PLAmi) = mn)\ Ri,^^ \'4,{n)) = Tr [ |;, m,)(/, m,|p«(n 



(17) 
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where, 



Ri,„n ^ ls<^\l,mi){l,mi\ (18) 

is a projection operator onto the eigenstates of Lz, and 

p(')(n)=Tr(^)[|V;(n))(V(n)|], (19) 

is the reduced state operator for the spin L at time n and Tr*-*-* denotes a trace over the factor space TCs- We have 
written out the exphcit expression (^) to emphasize that the probabihty of obtaining each mi value is associated 
with a projector onto a subspace of the factor space Tii. 

For reasons related to this fact (which we will make clear in later sections), we are also interested in examining the 
probability distributions associated with components of the total angular momentum J = S + L. The probability of 
obtaining a given mj value upon measurement of Jz is given by, 

PjA^j) = ^\{i^{n)\s,l,ms,mj - m^)]^, (20) 

iris 

where |s,Z,r7is,TOj — m^) is an element of the orthonormal basis (^). The probability Pj^inij) is associated with a 
projector onto a subspace of the full Hilbert space Ti.. The dimension of each subspace is given by the number of pairs 
{ms,mi) that yield a given value of nij = ms + mi. 

The classical probability distributions associated with dynamical variables are obtained by partial integration over 
the accessible phase space. In the case of Lz , the continuous marginal distribution is given by, 

P{Lz) = J J J dSzdcj),dcj)i p,{Sz.^s,Lz,(l)i), (21) 

where for notational convenience we have suppressed reference to the time-dependence. The marginal probability 
distribution for the total spin component Jz is obtained by integration subject to the constraint Sz + Lz — Jz, 

P{J.) = J J J j dSzd4>sdLzd4>i p,{Sz,(t>s,Lz,(j>i) 5{Sz+Lz - Jz). (22) 

These classical distributions are continuous, though their quantum counter-parts are intrinsically discrete. To con- 
struct a meaningful quantum-classical comparison it is useful to discretize the classical distrbutions by integrating 
the continuous probabilities over intervals of width Ti = 1 centered on the quantum eigenvalues. In the case of the 
component Lz, the quantum probability PlS^i) then associated with the classical probability of finding Lz in the 
interval [mi — 1/2, mi + 1/2]. This is given by 

PLirni)^ P{Lz). (23) 

Jmi-l/2 

Similarly, in the case of Jz, we compare each quantum Pj^{m.j) with the discrete classical probability, 

pmj+l/2 

Pl{m,)^ P{Jz). (24) 

Jmj-l/2 

In the following discussion of the numerical results we will emphasize that, for chaotic states, the steady-state 
shape of the quantum and classical distributions should be compared with the corresponding set derived from the 
microcanonical state. Our model is non-autonomous, but the spin magnitudes are conserved. The appropriate 
classical microcanonical measure is a constant on the accessible manifold V — x S^. This follows from the usual 
equilibrium hypothesis that all accessible microstates are equiprobable, where equiprobability is defined with respect 
to the invariant measure (H). This microcanonical density projected onto the L^-axis produces the discrete, fiat 
distribution, 

Pr:{mi) = {2l + l)-\ (25) 
However, projected along Jz, the microcanonical distribution is not fiat, but has a tent-shape. 
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pmc( ) ^ for \mA>l-s 

= 2[TT l"^^!^'-^- (26) 

In quantum mechanics, the equiprobabihty hypothesis implies that the appropriate microcanonical state is an equal- 
weight mixture. This microcanonical state, sometimes called a random state, is proportional to the identity in the 
full Hilbert space H = Hs ® Hi- It produces the same projected microcanonical distributions, i.e. ( p5| ) for and 
(p6|) for J2, as the classical microcanonical state. 



A. Mixed Regime Chaos 

We consider first a classical parameter regime (7 — 1.215, r = 1.1, and a = 5) for which the kinematically accessible 
phase space V is highly mixed. The chaotic region appears to be connected (all chaotic initial conditions have the 
same largest Lyapunov exponent \l = 0.04) and covers about half of the kinematic surface. A projection of only 
the chaotic initial conditions onto the plane spanned by Sz and Lz reveals large regular islands surrounding the 
stable parallel fixed points (±5*2, iLz), with chaotic regions spreading out from the unstable anti-parallel fixed points 
{±Sz,tLz)- a similar projection of the regular initial conditions shows points not only clustered about the parallel 
fixed points but also spread along the axis Sz = Lz- 

We now consider the time-evolution of quantum and classical states concentrated in the chaotic zone near one of 
the unstable anti-parallel fixed points, with initial centroids directed along 9{0) = (20°, 40°, 160°, 130°). The quantum 
dynamics are calculated using quantum numbers s = 140 and I — 154. As shown in Fig. ^ at early times both 
the quantum distribution PL^{mi) (solid line) and the corresponding classical distribution Pf^ (mi) (dots) remain 
well-localised. Their initial differences are not distinguishable on the scale of the figures. (The dots are shifted to the 
right by half of their width.) By time-step n = 20 both quantum and classical distributions have broadened to the 
system dimension and begin to exhibit noticeable differences. As shown in Fig. ||, around n — 100 the distributions 
have begun to settle close to an equilibrium shape. In Fig. || the successive time steps n — 199 and n — 200 show that, 
although both the quantum and classical distributions have relaxed very close to the same equilibrium distribution, 
the quantum distribution exhibits rapidly oscillating fluctuations about the classical steady-state. 

Both the quantum and classical equilibrium distributions (projected along Lz) show significant deviation from the 
microcanonical distribution (psf). This is also true of the distribution projected along L^, which has a different non- 
uniform equilibrium distribution than that observed when projecting onto Lz (see the left box of Fig. Uniform 
marginal distributions would be expected if the classical mapping was mixing, in which case arbitrary initial densities 
(with non-vanishing measure) would relax to the microcanonical distribution. Since the accessible kinematic surface 
has large KAM surfaces in this parameter regime, the coarse-grained classical equilibrium distributions are not ex- 
pected to be flat. An unexpected feature of the results is the observation that the shape of the equilibrium quantum 
distributions so accurately reflects the details of the KAM structure in the classical phase space. This feature is most 
striking in the case of the distributions projected along Jz (see the right box in Fig. ^). The steady-state quantum 
and classical probability distributions Pj^ (™j) and Pj^ {"^j) are both sharply peaked about rrij = 0. This equilibrium 
shape is much more sharply peaked than the tent-shape of the projected microcanonical distribution, PY'^{m,j), given 
by (|2^) and also plotted in the right box in Fig. ^. The important point is that the additional localization of the 
quantum distribution can be understood from a standard fixed-point analysis of the classical map |l3| : the presence of 
KAM surfaces arising due to the stability of the parallel fixed points prevents the chaotic classical spins from aligning 
in parallel along the z-axis. Most remarkably, we find that the steady-state quantum distributions accurately repro- 
duce this parameter-dependent structure of the mixed classical phase space even for much smaller quantum numbers. 
We examine how the accuracy of this correspondence scales with the quantum numbers in section [VI]|. 



B. Regime of Global Chaos 

If we hold a = 5 and r — 1.1 fixed and increase the coupling strength to the value 7 = 2.835, then all four of the 
fixed points mentioned above become unstable ||l^. Under these conditions less than 0.1% of the surface V is covered 
with regular islands; the remainder of the surface produces a connected chaotic zone with largest Lyapunov exponent 
Al — 0.45. We will sometimes refer to this parameter regime as one of global chaos since the kinematically accessible 
phase space is predominantly chaotic. 

The dynamics of the classical and quantum distributions are much simpler in this regime. We find that initially 
localised distributions, launched from arbitrary initial conditions, relax to the microcanonical distribution on a very 
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short time-scale. To demonstrate this, we consider the dynamics of an initial quantum state with s = 140 and I = 154, 
and a corresponding classical density, launched from 9{0) = (20°, 40°, 160°, 130°). Though the initial distributions are 
the same as in the mixed regime, by time-step n ~ 6 the quantum and classical distributions have already spread to 
the system dimension and begin to exhibit noticeable differences. By time-step n ~ 12 both distributions have relaxed 
very close to the microcanonical distributions. We plot the equilibrium quantum and classical projected distributions 
Pi^(mj) and Pj^{mj) in Fig. |^ for time-step n — 50. The projected classical distributions are early indistinguishable 
from the microcanonical forms, P£^f{mj) and P|"°(r7ij ), and the quantum distributions again exhibit small fluctuations 
about the classical distributions. We have found that these equilibrium quantum-classical differences asymptote to a 
non- vanishing minimum when the measure of KAM surfaces becomes negligible. These minimum quantum fluctuations 
reflect characteristic deviations from the microcanonical state that arise because the equilibrium quantum state is 
pure, whereas the microcanonical state corresponds to a random mixture. 



V. RATES OF RELAXATION TO EQUILIBRIUM 



In order to characterize the time-scale of relaxation to equilibrium it is convenient to study the time-dependence 
of a scalar measure that is sensitive to deviations from the equilibrium state. A conventional indicator of this rate of 
approach to equilibrium is the coarse-grained entropy, 

H = -J2Pd^P^■ (27) 

i 

Here the {Pi} stand for the quantum probabilities associated with projectors onto some basis of microstates {e.g. the 
projected distributions discussed in the previous section). The sum ( |27| ) is a standard measure of the information 
contained in a probability distribution and is sometimes called the Shannon entropy. 

The Shannon entropy has a number of useful properties. First, unlike the von Neumann entropy Tr[plnp], the 
Shannon entropy is basis-dependent. It reduces to the von Neumann entropy if the 'chosen' basis diagonalizes the 
state operator. However, this basis, or, more precisely, the set of projectors onto the (time-dependent) spectral de- 
composition of the state operator, does not necessarily correspond to a set of classically meaningful observables. Our 
main interest is to examine correspondence at the level of classical dynamical variables, so we consider probabilty 
distributions associated with projectors onto the eigenstates of classically well-defined operators. The classical coun- 
terparts to these probability distributions are associated with some fixed partioning of the phase space into cells of 
width h along the axes of the associated dynamical variable. 

Second, whereas the von Neuman entropy of the total system is constant in time (Tr[plnp] = since p =Up){ip\), 



the basis-dependent Shannon entropy may have time-dependence even if the quantum state is pure. Thus (27) may 
be applied to examine the rate of relaxation of either pure or mixed quantum states. It is in this sense that we use the 
term relaxation, although the time-evolution is unitary in the quantum model (and volume-preserving in the classical 
model). 

Given some fixed partioning of the phase space, if a classical state remains evenly spread through the phase space 
cells it occupies, and spreads through the phase-space exponentially with time, then an entropy like ( p7| ) should grow 
linearly with time. In this section we show that this argument holds approximately also for quantum states launched 
from a classically chaotic region of phase space. The actual rate of relaxation of the quantum states is accurately 
predicted by the classical entropy even for small quantum numbers. 

We demonstrate this behaviour by first considering the quantum entropy i?q[Jj] of the probabilities associated 
with the eigenvalues rrij of J^, i.e. the probabilities defined in (|2C|). The corresponding classical entropy, Hc[Jz], is 
calculated using the discrete classical probabilities (|2^). In Fig. ^~we compare the time-development of the quantum 
and classical entropies using quantum numbers s = 140 and I — 154. For these quantum numbers, the microcanonical 
{i.e. maximum) value of the entropy is H"^'^[Jz] = 6.2. In case (c), corresponding to a regular zone of the mixed 
regime {0{0) = (5°, 5°, 5°, 5°), 7 = 1.215), we actually see the greatest amount of difference between the quantum 
{Hq[Jz]) and classical {Hc[Jz]) entropies. Hq exhibits a quasi-periodic oscillation about its initial value whereas for 
He these oscillations eventually dampen. For smaller quantum numbers, and thus broader initial states. He dampens 
much more rapidly although Hq continues to exhibit a pronounced quasi-periodic behaviour. In case (b), with initial 
centroid 9{0) = (20°, 40°, 160°, 130°) set in a chaotic region of the equally mixed regime, both Hq and He oscillate 
about an initially increasing average before relaxing towards a constant value that lies well below the microcanonical 
maximum iJ'"°[Jz] = 6.2. This saturation away from the maximum is expected in the classical model since a large 
fraction of the kinematic surface is covered with regular islands and remains inaccessible. In case (a), corresponding to 
the regime of global chaos (7 = 2.835) and with the same initial state as (b), the entropies are nearly identical. Both 
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grow much more quickly than in the mixed regime case, roughly linearly, until saturating very near the maximum 
value. 

The quantum entropy is very well approximated by its classical counterpart also for smaller quantum numbers. In 
Fig. we display the growth rates of the quantum and classical entropies of the probabilty distributions associated 
with the observable for three sizes of quantum system (I — I — 22, I — 220) using the same parameters 
and initial condition as for data-set (a) in Fig. ^. In each case the quantum entropy is essentially identical to the 
corresponding classical entropy. The initial rate of growth is similar in each case, roughly linear, and of order the 
Lyapunov exponent, Al = 0.45. 

These results extend previous work demonstrating the that the widths of quantum states grow exponentially with 
time, on average, until saturation at the system dimension |^,^,^. Modulo the small quantum fluctuations, for 
both quantum and classical models we find that the subsequent relaxation to an equilibrium configuration occurs on 
the time-scale, 

trel-tsat+0(Xj}), (28) 

where t^at — A~^lnZ estimates the time it takes the initial coherent state to reach the system dimension. The 
exponent A^, is the exponent governing the growth rate of the state width H. The last term 0(A£^) approximates 
the additional time-required for the state to become more or less uniformly spread over the accessible phase space. 
In predominantly chaotic regimes we have found that A^ ~ Al, though in mixed regimes X^, is generally a few times 
larger than the largest Lyapunov exponent. 



VI. TIME-DEPENDENCE OF QUANTUM-CLASSICAL DIFFERENCES 

Before examing the scaling of quantum-classical differences with increasing quantum numbers, it is useful to deter- 
mine first their time-domain characteristics under the different types of classical behviour. Previous work has shown 
that quantum-classical differences for low-order moments, though initially small, grow exponentially with time when 
the classical motion is chaotic fl^Jl^ until the states approach the system size jl^. On this saturation time-scale 
those quantum-classical differences reach their maximum magnitude, but surprisingly this maximum was small, 0(K). 
More specifically, it did not scale with the quantum numbers. Of course, two distributions can be altogether different 
even when the differences between their means and variances are quite small, and therefore it is useful to examine the 
differences between the quantum and classical states in a more sensitive way. 

In this section we examine the time-dependence of bin- wise deviations between the quantum and classical probability 
distributions. For the observable L, this indicator takes the form. 



1 ' 

\ rni = -l 



(29) 



This standard deviation estimates typical quantum-classical differences on the scale of h along the L^-axis. Each 
interval is centered on a quantum eigenvalue. The P£ (m;) correspond to a measurement, or coarse-graining, of the 
classical density on an extremely fine-scale. 

In Fig. ^ we examine the time dependence of (t[Lz] for the same three classical sets of parameters and initial 
conditions displayed in Fig. ^. The initial value of (j[Lz] is generally not zero since it is not possible to match all 
the marginal distributions exactly in the case of the SU(2) coherent states The actual magnitude of the initial 
discrepancy depends on the angle between the axis of measurement, e.g. L^, and the direction of polarization of 
the initial state. For both chaotic states the differences initially decrease from their angle-dependent value and then 
increase until saturation at a steady-state value. This steady-state value is reached much later in the mixed regime 
(upper solid line), than in the global chaos regime (lower solid line). It occurs on the time-scale, trei, on which the 
underlying distributions have reached their steady-state configurations (modulo the quantum fluctuations). 

As shown in Fig. ||, the quantum-classical differences are actually largest for the regular state (dotted line) of 
the mixed regime (7 = 1.215) at both early and late times (relative to the relaxation time-scale). The steady- 
state magnitude of the differences for the global chaos regime (7 = 2.835) is significantly smaller than the typical 
magnitude for the mixed regime. However, for larger values of the classical perturbation strength 7, this average 
steady-state magnitude does not decrease further (with the quantum numbers held fixed) but has reached a non- 
vanishing minimum. The magnitude of the minimum steady-state fluctuations, Dl^ ~ 2 x 10~^, should be compared 
with a typical magnitude of the quantum and classical distributions, P^^ {mi) ~ 3 x 10~^. In the following section we 
examine how these fluctuations scale with increasing quantum numbers. 
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Above we have considered quantum-classical differences for observables (projectors onto subspaces) associated with 
the factor space 7i;. In this factor space the state is initially pure but becomes mixed as a result of dynamical 
interacions with the other subsystem. It is interesting to check if the dynamical behaviours of the differences are an 
artefact of this dynamical mixing. Therefore we consider also bin- wise quantum-classical differences for an observable 
(Jz = Sz + Lz) that acts non-trivially on the full Hilbert space Ti. ~ Hs <8i Hi. The quantum state in the full 
Hilbert space remains pure throughout the timc-cvolution. We construct the same standard deviation of the bin-wise 
differences between the quantum and classical probability distributions as above, 



= ^ [,(,^,)^,)] E[^^^K-)-PlK)F, (30) 

where rrij £ {l+s, l+s — 1, . . . ,~{l+s)}. In Fig. ^we compare a[Jz] in the same three classical regimes examined in Fig. 
^. Once again the regular state (dotted line) exhibits the largest quantum-classical differences, and the differences for 
both chaotic states (middle and lower solid line) grow to a steady-state value on the time-scale at which the underlying 
distributions relax to their equilibrium configurations. As above, the average value of the differences for 7 = 2.835 
(lower solid line) correspond to a non-vanishing minimum, that is, the average value does not noticeably decrease for 
larger values of 7. The minimum quantum fluctuations are again small when compared with the average height of 
the probabihty distribution, [2(s -I- /) -I- 1]"^ ~ 2 x 10"^. 

For 7 ~ 2.835, the measure of regular islands is already very close to zero and the classical system is nearly ergodic 
on v. Similarly, the quantum state is no longer constrained by any invariant classical structures but spreads almost 
evenly about the accessible Hilbert space. We find that the standard deviations of the quantum fluctuations that 
account for the equilibrium quantum-classical differences approach a non- vanishing minimum as the classical dynamics 
approach ergodicity on V. These equilibrium differences can not vanish (for fixed quantum numbers) because the 
total quantum state remains pure under the unitary dynamics, whereas the microcanonical equilibrium corresponds 
to an equal-weight mixture. 



VII. CORRESPONDENCE IN THE CLASSICAL LIMIT 



We now turn to an examination of the classical limit, /h —^ 00, where is characteristic system action. Since 
the quantum-classical differences grow to their largest values once the states have spread to the system size and 
subsequently fluctuate about this magnitude, we will examine the scaling of the differences in this late time-domain, 
that is, when the states have relaxed close to their equilibrum configurations. Moreover, these scaling results will then 
complement previous work that has focussed on correspondence at early times fl^ , in the Ehrenfest regime when the 
states are narrow relative to the system dimensions. 

We wish to determine if the standard deviation of the quantum-classical differences (defined in the previous sec- 
tion) decreases in magnitude with increasing quantum numbers. When comparing models with increasing quantum 
numbers, we hold the width of each probability bin fixed (at h = 1). Since the number of bins will increase with the 
quantum number, it follows that the height of the probability distribution in a given bin will also decrease. Conse- 
quently, we construct a scale- independent, or relative, measure of the bin-wise quantum fiuctuations by taking the 
ratio of the standard deviation to the average value of the probability distribution. For the observable Lz this takes 
the form, 

R[Lzin)] = £1^M = Ni a[Lz{n)]. (31) 

where the average value Pl^ = 1/(2/ + 1) = 1/Ni. If this relative measure approaches zero in the classical limit then 
the quantum probabilty distribution converges to the corresponding classical one in that limit. 

In Fig. ^ we consider typical equilibrium values of R[Lz{n)] plotted against We study the scaling using 

Ni because it is equal to the dimension of the factor space Hi and it is also proportional to the subsystem size 
Ni ~ 2|L|. We first consider a state launched in the global chaos regime (7 — 2.835, r — 1.1), with initial condition 
6{0) = (45°, 70°, 135°, 70°). The scatter of plus signs for each Ni — 21 + 1 value corresponds to time-steps n such that 
41 < n < 50. These time-step values are chosen because they occur well after the relaxation time trei — 6. In this 
regime the data exhibits very little scatter. A least-squares fit to the curve R = A/^/Nl + B yields a value for the 
intercept B that is consistent with zero [B = O.OOliO.OOl) and a slope of order unity {A — 1.032±0.02). An intercept 
consistent with zero implies that quantum-classical differences vanish in the classical limit, i.e. P]^^{mi) —^ P£ (w;) 
as I 00. This result is especially remarkable since we have considered the differences that arise given classical 
measurements which resolve the observable Lz with the rather extraordinary precision oi h — 1. 
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We next consider a state launched from the chaotic zone of the mixed regime (7 = 1.215, r = 1.1) with 9{Q) = 
(20°, 40°, 160°, 130°). The scatter of crosses in Fig. |l0| corresponds to time-steps 191 < n < 200 n, again chosen 
well after the relaxation time Uei for the range of quantum numbers considered. The scatter of quantum-classical 
differences at each Ni value is much more significant in this regime in which the equilibrium distributions reflect a 
much more complex phase space structure. However, the relative differences exhibit, on average, a similar dependence 
on the quantum numbers as in the predominantly chaotic regime. In this regime a least-squares fit to the curve R — 
A/Ni^/'^+B yields a slope of order unity {A = 3.39±0.15) but a negative value for the intercept {B = -0.017±0.009) 
within two-standard deviations of zero. A negative intercept is not physically meaningful (since i? is a positive definite 
quantity) and we assume it arises as a consequence of the statistical scatter in the data. Also plotted is the curve 

1 /2 

R = C/N^ , with slope C — 3.09 ± 0.04 also determined from a least-squares fit. Both fits are good, with reduced 
values of order unity. 

As we noted in the last section, the subsystem states do not remain pure, because of dynamically induced entangle- 
ment between the subsystems. Since the subsystem state ( p^ in the factor space is not pure, but highly mixed in 
the equilibrium time-domain, it is possible that the scaling with Ni that we observe is related to the purity-loss from 
this entanglement. Consequently, it is useful to examine the scaling of the quantum-classical differences for the total 
spin Jz- The operator Jz acts non-trivially in the full Hilbert space H. In this Hilbert space the system is described 
by a pure state vector at all times. In Fig. ^ we consider the scaling of the ratio, 

R[Jzin)] = ffiM = a[Jzin)], (32) 

where Pj^ = [2{s + I) + 1]-^ = Nj'^ is the average value of either distribution, versus the dimension Nj. Here Nj 
is the number of subspaces associated with distinct eigenvalues {rrij) of the quantum operator (Jz)- In contrast with 
iVj, Nj is not equal to the dimension of the corresponding Hilbert space, though it is a measure of the system size 
since Nj ~ 2|J|. The parameters and initial conditions shown in Fig. |ll] are the same as in Fig. |l^. The same fit 

1/2 

procedure as above, but applied to the function R = A/Nj + B, yields a value for B that is again consistent with 

zero {B = 0.00038 ± 0.0016) and a positive slope of order unity (A — 2.00 ± 0.04) in the predominantly chaotic 

regime (scatter of plus signs) . Thus the relative standard deviation for Jz also decreases as the square of the quantum 

numbers and fits to an intercept that is consistent with zero. This implies that the fluctuating quantum distributions 

approach the classical equilibrium, even for a few degree-of-freedom system, which is described at all times by a pure 

state. In a chaotic state of the mixed regime (scatter of crosses), the fluctuations are larger, and the same flt procedure 

as above gives {B = —0.016 ±0.012, A = 6.4 ±0.3), where the negative value for B lies within two standard deviations 

1/2 

of zero and is presumed to result from the statistical scatter of the data. Also plotted is the equation R = C /N-' , 
with C = 5.97 ± 0.06 determined from a least-squares fit. The fits to both equations are good, with reduced values 
of order unity. 



VIII. DISCUSSION 



We have shown that, in classically chaotic regimes, initially localised quantum states relax to an equilibrium config- 
uration that reflects the details of the classical phase-space structure. We find a remarkable degree of correspondence 
between the quantum and classical relaxation rates, even for small quantum numbers. Moreover, contrary to results 
obtained for the low-order moments [|l2|,^ , the degree of difference between the probability distributions is actually 
smaller for the chaotic states than the regular states. 

The equilibrium quantum distributions exhibit small rapidly oscillating fiuctuations about the coarse-grained clas- 
sical equilibrium. As the measure of regular islands on the classical manifold approaches zero, the quantum and 
classical equilibrium configurations approach their microcanonical forms, and the quantum fluctuations about the 
classical equilibrum approach a non-vanishing minimum. This minimum arises because we consider total quantum 
states that are pure, whereas the microcanonical configuration is produced by an equal-weight mixture. 

For the distributions associated with the subsystem observable L, the scale- independent standard deviation of these 
differences decreases as N^ ' where Ni = 2Z + 1 ~ 2|L| is the dimension of the factor space, and becomes vanishingly 
small in the limit of large quantum numbers (i.e. large spins). These results suggest that correspondence with classical 
Liouville mechanics emerges in the classical limit for time-scales much longer than the Ehrenfest time. 

A great deal of recent work has emphasized that the loss of purity resulting from interactions with a quantum 
environment removes characteristic quantum effects and improves the degree of quantum-classical correspondence 
J15|,p|p^. While this is certainly the case for small quantum systems, it has been further argued that these decoherence 
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effects must be taken into consideration to see the emergence of classical properties from quantum mechanics, even 
in the limit of large quantum numbers, if the classical motion is chaotic |5|l^. 

Since our model is comprised of interacting subsystems, initially separable pure states become entangled dynami- 
cally; the subsystem states (in each factor space) do not remain pure but become mixed. This entanglement process 
has an effect that is analogous to the process of decoherence. Hence one might suspect that the emergent classical 
behaviour that we have observed for the properties of the subsystem L may be strictly the result of a "decoherence" 
effect arising from entanglement with the other subsystem. To address this possibility, we have considered also the 
quantum-classical differences that arise in the probability distributions for a total system observable, Jz- In this 
case the quantum observables are projectors onto subspaces of the full Hilbert space, rather than merely a factor 
space. The quantum state in this full Hilbert space is not subject to any entanglement or decoherence and remains 
pure throughout the unitary time-evolution. We have found that the scale-independent standard deviations for these 
quantum-classical differences decrease as 1/ yJWj, where Nj — 2{s -|- ^) -|- 1 is a measure of the system size, Nj ~ 2|J|. 
where N = (2s -I- 1){21 + 1) is the dimension of the Hilbert space. The bin- wise quantum-classical differences become 
increasingly difficult to observe, in the limit of large quantum numbers, even for system observables that are isolated 
from the effects of decoherence. In this sense the process of decoherence is not necessary to produce quantum-classical 
correspondence in the classical limit. 
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FIG. 1. Quantum and classical probability distributions for Lz with / = 154 in chaotic zone of mixed regime 
(7 = 1.215, r = 1.1, a = 5). The dots are visible because they are shifted to the right by half of their width. The figure 
on the left is the initial state ( n = 0) and that on the right is at time-step n = 6. 
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FIG. 2. Same as Fig. |l| but for time-steps n = 99 on the left and n = 100 on the right. Both quantum and classical 
distributions have reached the system dimension and are relaxing towards equilibrium. 
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FIG. 3. Same as Fig. |l| but for n — 199 on the left and n — 200 on the right. The quantum distribution is fluctuating about 
a classical steady-state. 
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FIG. 4. Same as in the previous figure, but for PL^{mi) on tlie left and Pj^{mj) on the right, at time-step n = 200. Both 
Pjj,(mj) and Pj^{mj) are localised relative to the projected microcanonical distribution PJ^'^{mj). 
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FIG. 5. The equilibrium shapes of PL^imi) and Pj^{mj) at time-step n = 50 with I = 154 for a state launched in the 
global chaos regime (7 = 2.835, r = 1.1, a = 5). The quantum distributions exhibit small rapidly oscillating fluctuations about 
the projected microcanonical distributions. The classical distributions are not visible since the points lie within the fluctuating 
quantum data. 
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FIG. 6. Comparison of tlie quantum and classical entropies -ff [Jz] = — Pj^ {mj) log Pjz (f^j) for s = 140 and I = 154 in 
(a) regime of global chaos (7 = 2.835); (b) chaotic zone of the mixed regime (7 = 1.215); (c) regular zone of the mixed regime 
(7 = 1.215). 
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FIG. 7. Comparison of the quantum and classical subsystem entropies H[Lz 
system sizes in the global chaos regime of Fig. ^. 
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FIG. 8. Time dependence of the standard deviation (t[Lz] of quantum-classical differences ( |29| ) for states launched from 
a regular zone (dotted line) of the mixed regime (7 — 1.215), from a chaotic zone of the same mixed regime (middle solid 
line), and from the regime of global chaos (lower solid line, 7 — 2.835). The initial discrepancy is relatively large, but quickly 
decreases, and then increases until reaching an asymptotic equilibrium value. This occurs more slowly for the mixed regime 
case, for which the asymptotic value is also larger. In all cases s = 140 and I = 154. 
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FIG. 10. Scaling of relative quantum-classical differences in the equilibrum time-domain versus increasing system size. 
Scatter of crosses corresponds to time-steps 191 < n < 200, for a state launched in the chaotic zone of the mixed regime 
(7 = 1.215). Scatter of plus signs corresponds to time-steps 41 < n < 50, for a state launched in the global chaos regime. Data 
sets in both of these regimes are consistent with the scaling law R ~ Nc^^'^, where Ni = 21 + 1. 
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